Paramagnetic adsorbates on graphene: a charge transfer analysis 
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We introduce a modified version of the Hirshfeld charge analysis method and demonstrate its ac- 
curateness by calculating the charge transfer between the paramagnetic molecule NO 2 and graphene. 
The charge transfer between paramagnetic molecules and a graphene layer as calculated with ab 
initio methods can crucially depend on the size of the supercell used in the calculation. This has 
important consequences for adsorption studies involving paramagnetic molecules such as NO 2 ph- 
ysisorbed on graphene or on carbon nanotubes. 

PACS numbers: 68.43.-h, 73.20.Hb, 68.43.Bc, 81.05.Uw 



The use of carbon nanotube d 1 * 2 * 3 * 4 * 5 ^ and more 
recently graphene 7 as very sensitive gas sensors 
has stim ulated a lot of theoretical work on this 
subject l^ 1 11 ' 12 ' 13 1 14 1 1 ^ 1 Through ab initio calculations, 
one investigated the adsorption properties of these car- 
bon materials and looked for the mechanisms behind 
their good sensing capabilities. It appears that the key 
issue is the charge transfer from the gas molecules to 
the carbon surface. Such ab initio calculations were able 
to provide good qualitative agreement with experiment, 
e.g. whether the gas molecules act as electron donors 
or acceptors, but a large variation in the size of the 
doping was found between different theoretical calcula- 
tions. For example, the calculated charge transfer be- 
tween a (10,0) single- walled carbon nanotube (SWNT) 
and a NO2 molecule varies from 9 — 0.015e tcP^ — O.lOe, 
an order of magnitude difference! The size of the doping 
is, however, a crucial factor determining the sensitivity 
of the gas sensor. In this letter we demonstrate that im- 
portant reasons for the discrepancies are: i) the different 
sizes of the supercells used in these ab initio calculations, 
and ii) the theoretical charge analysis method that has 
been used. We will show that the former is a very deci- 
sive parameter when the adsorbing molecules are param- 
agnetic. First we examine critically different theoretical 
charge analysis methods. Then we will perform ab initio 
calculations with different supercells to investigate the 
charge transfer dependence on the size of the supercell. 

All our ab initio calculations were done within the 
density funcional theory (DFT) formalism using the 
ABINIT^ software package within the local spin den- 
sity approximation (LSD) and with Troullier-Martins 
pseudopotentialsP^ We used a plane wave basis set with 
an energy cutoff of 816 eV, which was tested to give con- 
verged results for all the properties studied in this letter. 
Different graphene supercells ranging from 2 x 2 to 6 x 6 
were implemented. For the sampling of the Brillouin zone 
(BZ) we used Monckhorts-Pack (MP) grids for the dif- 
ferent supercells equivalent to a range from 12 x 12 x 1 
to 48 x 48 x 1 points for a single unit cell. 

To calculate the charge transfer from a molecule to 
a surface, one needs a physically meaningful and trans- 
parent approach to divide the electron density between 
them. A variety of methods has been developed for this 



purpose, such as the Mulliken's 18 or Bader's 19 atoms in 
molecules approaches. We introduce a method based on 
the Hirshfeld 20 method and in particular the recently 
extended version of this approach, the Hirshfeld-P^ 
method. The simplicity and low cost of computation time 
makes the Hirshfeld method a very powerfull approach in 
DFT calculations, and we will demonstrate that its ex- 
tended version is more accurate than the Bader method 
and other Hirshfeld-based methods. In the Hirshfeld 
charge analysis, the total electron density p is divided 
between the different atoms of a system, according to the 
density of the neutral atoms, p A , in free space which build 
up this system: Q A = J (p° A (r)/Y,A' Pa' ( r ))p( r )^ r > with 
Qa the charge on atom A (in units of e). The Hirshfeld- 
I method does not use the density of neutral atoms in 
free space, but uses the density of charged atoms in- 
stead. These charges on the different atoms are deter- 
mined through an iterative process. One starts with the 
simple Hirshfeld method and determines the charge per 
atom. Then this charge is placed on the different atoms 
in free space and their densities are used to divide the 
density again among the atoms as with the simple Hirsh- 
feld method. This procedure is repeated until the charges 
on the atoms are converged. In our approach we allow 
our system to relax completely and calculate the total 
density. Then we calculate the density of the molecule 
and the graphene layer separately in the same configu- 
ration as in the total relaxed system. These densities 
can now be used in the Hirshfeld-I method instead of the 
separate atomic densities. Such an approach gives charge 
transfers that are more physically meaningful as we will 
demonstrate in the following. 

To test this charge analysis procedure we investigate 
the adsorption of the paramagnetic NO2 molecule (with 
magnetic moment M = 1^) on a 4 x 4 graphene su- 
percell. We calculate the charge transfer of the NO2- 
graphene system using different theoretical approaches: 
we first use the simple Hirshfeld method, then we change 
the atomic densities into a molecular and a graphene den- 
sity and use the same procedure; finally we put charges on 
NO2 and graphene in order to use them in the Hirshfeld-I 
charge analysis. The Hirshfeld methods that use molec- 
ular densities instead of atomic densities will be referred 
to as the modified Hirshfeld and modified Hirshfeld-I 
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method. For comparison, we also calculated the charge 
transfer based on the Bader charge analysis. The charge 
transfer results are given in table [TJ 

TABLE I: The charge transfer from graphene to NO2 calcu- 
lated with different methods. 

method charge transfer (e) 
Hirshfeld -0.099 
modified Hirshfeld -0.161 
modified Hirshfeld-I -0.181 
Bader -0.212 
magnetic moment -0.182 



The advantage of using a paramagnetic molecule to 
demonstrate the accurateness of this charge transfer anal- 
ysis method is that we can have a good estimate of the 
charge on the molecule through some physically relevant 
properties of the system, like the magnetic moment: ex- 
tra charge in the partially occupied molecular orbital 
(POMO) of NO2 will lower the magnetic moment of the 
NO2 molecule. Extra charge on the graphene layer, on 
the other hand, will not effect its magnetic moment be- 
cause pure graphene is diamagnetic 22 . This implies that 
we can simply use the lowering of the magnetic moment 
of the total system (i.e. the difference from to esti- 
mate the charge transfer between the NO2 molecule and 
graphene. This is very accurate if there is not too much 
hybridization between the molecular orbitals of NO2 and 
the graphene orbitals. This is i ndeed the case for NO2 
since it physisorbs on graphene! 13 * 14 ! The value of the 
charge transfer extracted from the lowering of the mag- 
netic moment is — 0.182e (see table |l|, which is very 
close to the charge transfer obtained from the modified 
Hirshfeld-I method. The use of molecular densities in- 
stead of atomic ones leads already to a significant im- 
provement of the simple Hirshfeld method and the it- 
erative procedure makes the calculated charge transfer 
almost equal to the one extracted from the change in 
magnetic moment of the total system. The small differ- 
ence of the order of O.OOle is caused by a small charge 
transfer due to orbital hybridization. This causes no 
change in the magnetic moment, but it is noticeable in 
the modified Hirshfeld-I charge analysis. Note also that 
the charge transfer from the Bader analysis does not cor- 
respond with the charge transfer from the lowering of 
the magnetic moment (see table |l| . Thus the modified 
Hirshfeld-I method gives a good approximation to the 
charge transfer in the adsorption process of paramagnetic 
molecules on graphene and, if the molecules are not para- 
magnetic, it is probably also more accurate to use this 
modified Hirshfeld-I charge analysis instead of e.g the 
Bader charge analysis for a trustworthy determination of 
the charge transfer. 

Another element of support for the modified Hirshfeld- 
I method is its geometrical interpretation. When we 
put an extra charge on a free standing paramagnetic 
molecule, it will be placed at the POMO. In the case of 



a NO2 molecule, this is an anti-bonding orbital which 
means that the bond length between the N and O atoms 
becomes larger when the POMO becomes filled. This 
fact can be used to estimate the charge transfer of a 
physisorbed NO2 molecule to graphene by comparing 
the bond length of the relaxed physisorbed molecule 
with a charged free standing NO2 molecule (see Fig. [I]). 
From Fig. [T] we extract a charge transfer of — 0.206e, 
which is close to the charge transfer obtained by the 
modified Hirshfeld-I method (and not close to the simple 
and modified Hirshfeld methods). 




bond length a (A) 

FIG. 1: (Color online) The dependence of the bond length of 
a NO 2 molecule on the extra charge put on the molecule. 

Next we will discuss how the charge transfer for the 
case of paramagnetic molecules depends on the size of 
the supercell used in the ab initio calculations. The sim- 
ulation of adsorption processes at surfaces is necessarily 
restricted to finite supercells. In most cases this gives a 
reasonably good approximation for e.g. charge transfer 
analyses in physisorption processes, but in some cases, 
when e.g. the adsorbates are paramagnetic, a different 
kind of charge transfer mechanism dominates the charge 
transfeJ^^T leading to very different results. We will 
use the NO2 molecule again to show how sensitive the 
charge transfer can be with respect to the size of the 
used supercell. We make use of 4 different supercells, 
2 x 2,3 x 3,4 x 4, and 6x6 and several MP-grids for 
sampling the BZ. Fig. [2] shows that the charge transfers 
are converged reasonably well for a 24 x 24 x 1 MP-grid. 

There is a pronounced dependence of the charge trans- 
fer on the number of atoms in the simulated graphene 
layei^ which can be explained by lo oking at the den- 
sity of states (DOS) of the total systenP™ ( see Fig. [|]). 
The lowest unoccupied orbital (LUMO) of NO2 is be- 
low the Dirac-point 24 of graphene, which causes a charge 
transfer to the molecule. This charge transfer depends 
clearly on the number of electronic states between the 
Dirac-point and the LUMO of the NO2 molecule, which 
depends linearly on the number of carbon atoms in the 
supercell. Therefore, we may expect that if one takes a 
graphene supercell that is large enough, one would even- 
tually get a charge transfer of one electron!^ However, it 
is clear from Fig. [2] that there is no simple linear depen- 
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FIG. 2: (Color online) Supercell dependence of the charge 
transfer between a NO2 molecule and graphene for different 
MP-grids. The inset shows a 4 x 4 graphene supercell. 



dence. The position of the LUMO in the DOS depends 
on the filling of this orbital: the more electrons in the 
orbital the more difficult it gets to put another one in it 
and this translates in a shift of the LUMO towards the 
Dirac-point. In Fig. [4] we show this shift as a function of 
the number of carbon atoms in the supercell. It is not 
obvious that the orbital becomes totally filled (a trans- 
fer of one electron per paramagnetic molecule) before it 
coincides with the Dirac point, where the shift will stop 
and the charge transfer will be converged at a value of 
less than one electron per molecule. 

Thus one can only obtain a quantitative meaningful 
value of the charge transfer between a paramagnetic 
molecule and a graphene layer in two cases: i) when the 
supercell used is large enough to get a charge transfer 
of one electron per molecule, or ii) the LUMO of the 
molecule coincides with the Dirac-point and the charge 
transfer has converged. In the case of the NO2 molecule, 
one needs a supercell that is much larger than the 6x6 
supercell we were able to use due to computational 
limitations. However, from Fig. [4] we find that if the 
LUMO stays partially full, it takes more than thousand 
carbon atoms to let it coincide with the Dirac-point. On 
the other hand, from Fig. [2] one may deduce that it is 
not unlikely that a charge transfer of le is reached before 
the supercell contains 1000 carbon atoms. From this we 
can conclude that probably the LUMO of NO2 will be 
completely filled before it reaches the Dirac-point. This 
is compatible with the experimental observation 7 that 
the charge transfer equals one electron charge per NO2 
molecule. 

To conclude, we examined different methods to calcu- 
late the charge transfer between a paramagnetic molecule 
and a graphene layer and found that a modified Hirshfeld- 
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FIG. 3: (Color online) Spin-polarized DOS for NO2 adsorbed 
on graphene calculated with a 60 x 60 x 1 MP-grid for a 4 x 4 
supercell. The LUMO of the NO2 molecule and the Fermi- 
energy are indicated. The inset shows the orbit als of the 
NO2 molecule (the LUMO) and graphene between which the 
electron transfer takes place. 
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FIG. 4: (Color online) Distance, AE, of the N0 2 LUMO to 
the Dirac-point with respect to the number of carbon atoms 
used in the supercell. The inset shows a close-up of the total 
DOS of graphene with an adsorbed NO2 molecule around the 
Dirac-point. 



I method gives the most accurate and physically mean- 
ingful results. We used this method to investigate the 
dependence of the charge transfer on the size of the super- 
cell in ab initio calculations and found that charge trans- 
fers involving paramagnetic molecules are, in a nontrivial 
way, very sensitive to the supercell size. We applied our 
results to the adsorption of NO2 on graphene and showed 
that our ab initio calculations are compatible with the 
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experimentally found charge transfer of leJ In contrast 
to the claim in Ref. [T3J we found that a charge transfer 
of le is not always realized for adsorption of a param- 
agnetic molecules on graphene. The present results are 
also valid for single wall carbon nanotubes because their 
adsorption properties are similar to those of grapheneP^ 
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